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The energy associated with bosonic and fermionic pairs of topological spin defects in doped an- 
tiferromagnetic quantum spin- 1/2 square lattice is estimated within a resonating valence bond sce- 
nario, as described by a i-i'-J-like model Hamiltonian, plus a t±, responsible of a three-dimensional 
screening of the electrostatic repulsion within the bosonic pairs. For parameters appropriate for 
monolayered high-Tc superconductors, both fermionic and bosonic pairs show x 2 — y 2 symmetry. 
We find a critical value of doping such that the energy of the bosonic pairs goes below twice the 
energy of two fermionic pairs at their Fermi level. This finding could be related to the onset of 
high-r c superconductivity. 
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I. INTRODUCTION 

Since the discovery of high-T c superconductiv- 
ity (HTSC) in La2CuC>4?i a vast amount of work has 
been done on slightly doped quasi-two dimensional (2D) 
antifcrromagnetsi 2 ! 3 ! 4 ! 5 ! 6 It has been observed that these 
materials display very unusual properties, with a rich va- 
riety of temperature-doping phases diagram. Specifically, 
away from the overdoped side, the cuprates do not ap- 
pear to be a Landau Fermi liquid. For instance, they 
should be considered as doped Mott insulators. However, 
the theoretical status of the field has been largely phe- 
nomenological and controversial^ As far as we know, 
there is no consensus on the origin of the superconduc- 
tivity nor on the pseudogap phase. Therefore, finding a 
microscopic mechanism for HTSC still is an open prob- 
lem. 

Local-density approximation (LDA) and generalized 
gradient approximation (GGA)£ to density functional 
theory have been used so far to rationalize the electronic 
structure of HTSCs. Although LDA is a useful tech- 
nique for some materials, it has been shown that both 
LDA or GGA are not appropriate for antiferromagnetic 
materials because they tend to yield a metallic ground 
state with incorrect delocalized spin density and band 
ordering^**! This is attributed to an extreme nonanalytic 
and nonlocal behavior of density functional theory as the 
particle number is changed^Siii implying the need for 
a self-energy correction, or at least an orbital depen- 
dent potential to obtain a realistic description of band 
gaps. To overcome such a problem, different semiempiri- 
cal corrections to LDA hav e been proposed so far as, i.e., 
LDA+SIC (Refs. ll2ll3ll4l and LDA+l/. 15 : 16 : 17 

An alternative approach to the electronic structure of 
the HTSCs is based on the use of model Hamiltonians 
that aim to incorporate the essential physics into a few 
parameters. It is generally accepted that electron cor- 
relation is important for HTSC. Furthermore, it is well 



known that the (covalent-structure) valence-bond (VB) 
model or, equivalently, the Heiscnberg Hamiltonian in- 
cludes most of the electron correlation. Thence, early 
in 1987, Andersoni^ proposed that the important fea- 
tures of the undoped HTSC parent compounds can be 
described by a Heisenberg Hamiltonian on a two-dimen- 
sional square lattice with one electron per site. Mean- 
while, Emerjii2i2£ proposed a three-band Hubbard model. 
Unfortunately, the number of parameters of a three-band 
Hubbard model turns to be too large. Therefore, Zhang 
and Rice 21 - proposed a simplification of the three-band 
Hubbard model into the well known t-J, which implicitly 
includes the 0(p)-Cu(d) hybridization and recovers the 
initial effective one-band description of Anderson. Since 
an appropriate parametrization is essential for the pre- 
dictive capability of model Hamiltonians, much progress 
has been achieved on the high-level ab initio computation 
of reliable appropriate parameters using only the crystal 
structure as external input £2/22/21 

Stimulated by Anderson's suggestion^ a renewed in- 
terest of low-dimensional quantum spin-1/2 antiferro- 
magnetic systems emerged. According to the Lieb and 
Mattis theorem 25 the ground state for the undoped half- 
filled bipartite system must be a singlet. Therefore, 
the appropriate ground-state wave function could have 
a resonating-valence-bond (RVB) character. It was soon 
pointed out that short-range RVB wave functions ex- 
hibit topological long-range order 2£/2L22£2\ Furthermore, 
recently 30 topological order for superconductors has also 
been claimed, away from truly microscopic models, mak- 
ing use of bosonic theories of the quantum Ginzburg- 
Landau form. In Ref. 28 Klein and collaborators investi- 
gated the short-range RVB wave functions within a dimcr 
coverings approximation for the square lattice and found 
that the dimer-coverings show a type of long-range spin- 
pairing order (LRSPO). Using arguments based on the 
LRSPO they predicted a per-site energy e oc <5 2 , where 
8 is the deviation of the local LRSPO with respect to 
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the LRSPO of the ground state. Furthermore, topolog- 
ical spin defects (TSDs), namely a site that is not spin 
paired to a singlet, or a hole in hole-doped superconduc- 
tors, or a doubly occupied site in electron-doped mate- 
rials, were assimilated to Bloch walls separating phases 
with a difference in LRSPO of ±1. It was argued that, at 
a longitudinal distance ~ p past the TSD on the less sta- 
bilized side, the defect should also presumably have only 
spread out a transverse distance ~ p, so that 5 ~ 1/p, 
and e oc 1/p 2 . Therefore, the energy contribution from 
all the sites of the given longitudinal distance past the 
TSD is pAe ~ 1/p. Summation over all the sites up to 
a given distance thence gives an energy cost of ~ In p. 
When the TSD are charged, it was also suggested that 
this long flat attraction ~ In p along with the screened 
repulsion exp{— ap} / p could lead to a weakly bound pair. 
Recently^ a linear relationship between T c and the J/t 
ratio, as obtained from high-level ab initio calculations, 
was found. It was argued that such a linear relation arises 
from the LRSPO mechanism previously suggested^ Fur- 
thermore, the so-called t-J Hamiltonian for the cuprates 
seems to be pointing to the right direction. 

The existence of a LRSPO for more general RVB wave 
functions has been proven for ladderlike quantum spin- 
1/2 antiferromagnetic systemsi2Li2£ Most of the consider- 
ations associated to the existence of this LRSPO for the 
ladderlike quantum spin-1/2 antiferromagnetic systems 
are readily applicable to the square lattice. In particu- 
lar, bound pairs of TSD are predicted to occur. How- 
ever, as far as we know, the energy of such a pair of 
TSD as a function of the distance has not been obtained 
yet. Even more, arguments based on the LRSPO alone 
cannot decide if vacancies (doubly occupied sites), let's 
say charge-wearing TSDs, organize themselves as bound 
pairs of two charge-wearing TSDs, as bosonic-character 
pairs, or each charge-wearing TSD would bind to a non 
spin-paired spin, let us say a spin-wearing TSD, leading 
to a fermionic-character pair. 

Here focus is on the energy associated with these 
bosonic and fermionic pairs as described by symmetry- 
adapted extended wave functions. We find that the 
fermionic pairs are favored for low doping levels, but the 
Fermi level increases with doping while the energy of the 
bosonic pairs lowers. At a critical doping the energy of 
the bosonic pair goes below the energy of two fermionic 
pairs at the Fermi level, suggesting the pairing of charge- 
wearing TSDs, and hence providing a microscopic mech- 
anism for HTSC. 

The description of these bosonic and fermionic pairs 
is based on a t-t'- J- like model Hamiltonian H = Hi + 
Hj+Ht+Hf, where Hi is the energy associated with the 
ionization potential for hole-doped materials, or the en- 
ergy associated with the electron affinity for the electron- 
doped systems. The Hj is the well known nearest- 
neighbor Heisenberg Hamiltonian, 

Hj = J Sr-Sr/, (1) 

<R,R') 



where Sr is the spin operator for the spin on the site 
R, and (R, R') means that R and R' are nearest neigh- 
bors. The nearest- and next-nearest-neighbor hopping 
contributions to the Hamiltonian are, respectively, 

H t = -t ( c L c R'tf + C R CT C ^) 

(R,R'> <r 

x (1 - h Rs ) (1 - hn>a) , (2) 

H t < = *'<(R,R')> 51 ( C R-^ CR ' CT + c R'^ CRct ) 

((R,R'» 

x (1 - n Rff ) (1 - riRff) , (3) 

where c Rct (cr ct ) creates (destroys) an electron on site R 
with spin a — a, f3. The double occupancy is avoided by 
the factors 1 — riRg- , where n-Ra is the number operator on 
site R with spin a = /3, a. The summation on ((R, R')) 
means that R and R' are next-nearest neighbors. The 
hopping integral i'^ R depends on the number of holes 
within the plaquette. Such a model is known to repro- 
duce the low-energy spectrum of the three-band Hubbard 
models Here, we use the parameters obtained by high- 
level ab initio calculations using only the crystal struc- 
ture as external input^S We approximate the screened 
electrostatic repulsion within the bosonic pair by the 
Yukawa potential^i^ the screening agent being the gas 
of the fermionic pairs. The three-dimensional (3D] char- 
acter of the screening is taken into account by an inter- 
layer hopping integral t±. The Heisenberg part of the 
energy associated with a pair of static TSDs is estimated 
by the dimer-covering- counting approximation2£i2L2£ on 
w x L antiferromagnetic quantum spin-1/2 square lat- 
tice, with w = 4, 6, ... , 20, L — > oo, and cyclic boundary 
conditions in both directions. Counting of the dimer- 
covering configurations has been achieved by a transfer- 
matrix techniquei2ii22 

This paper is organized as follows: In Sec. [H] we re- 
view the main concepts about LRSPO, 31 -jjH describing 
the scenario where the TSDs are located. In Sec. IIIH thc 
energy per site of the half-filled ground state, and the 
gain of the Heisenberg energy associated with a static 
pair of TSDs is estimated. In Sec. IIVI the symmetry- 
adapted extended wave- functions appropriate for bosonic 
and fermionic pairs moving in a Cu02 layer will be de- 
fined, and the energy bands will be obtained. From the 
two-fluids equilibrium condition, the critical doping p c 
for the onset of pairing to bosonic pairs among charge- 
wearing TSDs is obtained. Finally a summary and the 
conclusions are given in Sec. [3 

II. LONG-RANGE SPIN-PAIRING ORDER 
AND TOPOLOGICAL SPIN DEFECTS 

From the Lieb and Mattis theorem^ it is well known 
that for bipartite spin systems a maximally-spin-paired 
ground state is expected. In particular, at half filling, for 
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FIG. 1: A fragment of a VB configuration for a w — 8 square 
lattice. Each arrow represents a spin-pairing (SP) to a singlet 
between a spin on a site in the sublattice A (circles) with a 
spin on a site in the sublattice B (stars). Below each bound- 
ary, i.e. the dashed lines running parallel to the rungs, it ap- 
pears the net count of arrows, D, penetrating this boundary. 

ladder-like systems, with equal number of sites in the A 
and B sublattices, the ground state is a singlet. Singlet 
states can be achieved by configuration interaction (CI) 
among covalent VB configurations or RVB. For instance, 
a linearly independent set of VB configurations can be 
achieved by pairing to a singlet each spin in the sublattice 
A to a spin in the sublattice B, not necessarily one of its 
nearest neighbors (see Fig. . 



It is known^i^ that any (covalent) VB configuration 
exhibits a LRSPO related to the local (at boundary) 
array of SPs penetrating any boundary /„ (see, for in- 
stance, Fig. 0. The parameter associated with the LR- 
SPO, D, can take w + 1 different relevant values. The 
shape of the boundary selected to define the different 
w+ 1 values of D is quite arbitrary, though when w=even 
and the boundary is chosen to run parallel to rungs, the 
w + 1 different values of D are: 

D = 0,±1,...,±|. (4) 

This LRSPO allows to separate the set of VB config- 
urations in different subsets. Since two singlets from 
different subsets must be different repeatedly at every 
position along the ladder, they are asymptotically or- 
thogonal and non interacting via any interaction medi- 
ated by a few-particle operator. Then the matrix of the 
Hamiltonian asymptotically block-diagonalizes, so con- 
figurations belonging to different subsets do not mix in 
the CI sense. Thus D may be taken as a long-range or- 
der parameter labelling the eigenstates of the D block. 
Under low-frustration conditions, the expected ordering 
of the lowest-lying energy Ed from the different blocks 
is 

E < Ex < ■■■ < E w/2 , (5) 



with Ed = E-d- 

Now, half-filled excited states or slightly doped states 
are analyzed via TSDs. There are different types of exci- 
tations conceivable from a maximally- spin-paired ground 
state. Let us say, preserving half filling (one electron per 
site), there are primarily spin excitations. In this case, 
two spin-wearing TSDs, one in the sublattice A and the 
other in the sublattice B, are obtained by breaking one 
SP to form a triplet state. Away from half-filling, there 
are low-energy spin and charge excitations. For instance, 
removing (adding) one electron produces two sites that 
cannot be SP, a charge-wearing TSD and a spin-wearing 
TSD, one in the sublattice A and the other in the sub- 
lattice B, the ladder becoming a doublet. In this case 
hopping terms must be retained in the Hamiltonian and 
the so-called t-J model or different extensions that incor- 
porate either next-nearest-neighbor hopping t' or electro- 
static repulsion have been employed so far. Thence, the 
doublet is a weighted superposition of VB configurations 
with a spin- wearing TSD and a charge- wearing TSD lying 
in different sublattices. Still, going up in the hierarchy 
of Hamiltonians, the Hubbard or even a more general 
Hamiltonian must be considered. In this case, still an- 
other type of excitations (though presumably of higher 
energy if a Heisenberg-type Hamiltonian is assumed to 
govern the lowest-lying region of the spectrum) can be 
produced relaxing the single-occupancy constrain. This 
leads to the ionic states, i.e., states with at least a pair of 
sites, one doubly occupied and the other empty, namely 
one negatively charge-wearing TSD and one positively 
charge- wearing TSD. 

Of special relevance here is how the LRSPO is dis- 
rupted by a TSD (see Fig. 0). For instance, a TSD in 
a site [n,i], n indicating the rung and i the leg, can be 
seen as a domain wall on the rung n which separates the 
lattice in two sectors with associated left, Di, and right, 
D r , order parameters. When we choose the sublattice A 
as formed by the set of sites [to, j] with to + j=even, 

A = A - (-1)" +I . (6) 

Thence, to fulfill boundary conditions TSDs must appear 
by pairs, one TSD in the sublattice A and the other in 
the sublattice B, to ensure AD = from the left to the 
right of the pair. Such a pair define an intervening re- 
gion with AD — ±1 with respect to the LRSPO D of 
the host (see Fig. [2J. Then, away from half-filling, it 
may be conceivable an intervening region limited by two 
charge-wearing TSDs, or a charge-wearing TSD and a 
spin-wearing TSD (provided that the doping is not so 
strong as to preclude a maximally-spin-paired ground 
state). In particular, when placing a pair of TSDs above 
the ground state (D=0), the order parameter of the in- 
tervening region will be |Z? P | = 1, which from Eq. JSJ) is 
expected to have higher associated energy per site. This 
indicates that the pair of TSDs should try to remain as 
close as possible. Thus, bound pairs of TSDs are pre- 
dicted to occur. To show that this is the case is one of 
the concerns of the present paper. 
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FIG. 2: A fragment of a VB configuration for a w — 8 square 
lattice containing a pair of TSDs, one in the sublattice A 
(white circle), and the other in the sublattice B (white star). 
Notice that this VB configuration shows LRSPO D = ev- 
erywhere but in the intervening region defined by the pair of 
TSDs, with D = 1. 



III. HEISENBERG ENERGY OF A STATIC 
PAIR OF TSDS 

Within the dimcr-covering- counting approximation 
the resonance energy, E r {w,D) in units of J, i.e., the 
energy correction below the energy of a single dimer- 
covering structure, depends on the configuration inter- 
action among the different dimer-covering configurations 
with LRSPO D. When an equally weighted wave func- 
tion is considered, it has been arguec&SLiS that one 
might consider this energy lowering to depend solely on 
the dimension of the space spanned by the appropriate 
dimer-covering configurations. Let Nd(w) be the num- 
ber of linearly independent dimer-covering configurations 
with the LRSPO D. Since Md{u)) is multiplicative in 
terms of a break up into subsystems while the energy is 
additive, such a functional dependence should be of the 
form 

E T {w,D) m -ChxN D {w), (7) 

The values Md (w) can be easily obtained by a transfer- 
matrix technique^ Let us start computing Nd{w) for a 
maximally spin-paired half-filled system. Let us analyze 
from a local point of view the dimer-covering singlets. We 
can identify the dimer-covering local states according to 
which legs have a pairing across the f n boundary. In the 
present case it can be seen that, for each boundary, there 
are 2 W different local states, |e n j) (/ ranging), which can 
be classified according to the value of D, \e®j). 

Proceeding from the left to the right, from the bound- 
ar y fn-i to /„, a dimer-covering- counting matrix, /C„, 
is defined as (e n -u\lC n \e n j) = the number of different 
ways \e n j) can succeed |e n _ij). Then, the number of 



dimer-covering states in a I? subspace is 

A^( W )=^(e^|/C 1 /C 2 ---/C L |e? 7 ). (8) 

Since D n — D n+ \ for any dimer-covering singlet, K, n is 
a block-diagonal symmetric matrix that does not depend 
on n and we can omit this subindex. For L — > oo, the 
highest eigenvalue A w d of the D block JCd dominates, 
and 

N D (w)^A L wD . (9) 

Therefore, 

E D (w) ps wLe + E r (w, D) w wLe — CL lnA^rj, (10) 

where Eq = —0.375 is the energy per spin of a single 
dimer-covering configuration. Sinc o 31 i 32 

A wD >A wD , when \D\ < \D'\, (11) 

with K w o — A^ira, the Heisenberg energy for the half- 
filled ground state belongs to the subspace D=0, as sug- 
gested by Eq. J3J), and can be approximated (in units of 
J) by 

E (w) w wLe - CL\nA w0 . (12) 

C is a fitting parameter independent of the structure 
to some degree. The value of C for the nearest-neighbor 
isotropic Heisenberg model has been determined for a 
class of benzenoid hydrocarbon^ (with C=0.5667) and 
for finite square-lattice fragments 3 ^, (with C=0.75), by 
fitting the logarithm of the dimer-covering count to the 
resonance energy calculated from an equally weighted 
dimer-covering wave function. For the spin-1/2 square 
lattice, more general RVB approximations suggest a 
rough estimate ofC = 0.94 ±0.19. Here C is fixed 
to yield a reasonably good estimate of the ground-state 
Heisenberg energy of the half-filled square lattice. Ta- 
ble [I] summarizes the ground-state energy per site for 
w = 4, 6, . . . , 20 and its extrapolation to ro — * oo. We 
use C = 1 from here on, since C w 1.0083 yields the 
ground-state energy of Liang et alm& 

When adding a TSD to a Cu02 layer the transfer ma- 
trix K, across the defect must be substituted by the ap- 
propriate /Cr, where R is the vector position of the TSD. 
Therefore, the number of dimer-covering configurations 
when adding a pair of TSDs to the half-filled ground 
state, located, respectively, at and [m,j], with non- 
negative m and j, with m + j=odd, is 

M[m,j]{w) = (A w0 |/C /C m - 1 /C [mj] /C L - m - 1 |A tu0 ) 

w A^ m " 1 (A w0 |/C /C™- 1 /C [mj] |A to0 ). (13) 

Thence, the Heisenberg energy (in units of J) associated 
with a pair of static TSDs separated [m, j], m + j=odd, 
with respect to the energy of the half-filled ground state 
is 

e [mJi (w) « -2e + In "° (14) 
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TABLE I: The ground-state resonance energy per site in units 
of CJ (— ]hA w q/w), and the extrapolation to w — *• oo. In the 
third column the ground-state energy per site in units of J as 
obtained when taking C = 1 . 
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IV. TWO-DIMENSIONAL EXTENDED WAVE 
FUNCTIONS 



The bosonic or the fermionic pairs are far from be- 
ing static. The hopping terms of the Hamiltonian allow 
any charge- wearing TSD to move while the exchange part 
mixes up all the VB configurations. Therefore, the appro- 
priate wave function must be a weighted superposition of 
all possible static configurations, fulfilling translational 
and point group symmetry conditions. Thence, the wave 
functions for both bosonic and fermionic extended pairs 
of TSDs should be invariant under the operations of the 
factor group isomorphic to the C/± v (4 mm) group, as ob- 
tained by factorizing the full group into the translational 
subgroup and the planar subgroup. Thence, there can be 
conceivable extended wave functions with symmetry S, 




Ai, totally symmetric, (x 2 + y 2 ), 

A2, antisymmetric under the four reflections, 

B±, antisymmetric under Cf, cr x ± y , (x 2 — y 2 ), 

B2, antisymmetric under C 4 , a x , cr y , (xy). 



(15) 

Therefore, symmetry-adapted extended wave functions 
for both the fermionic pairs and bosonic pairs can be 
written as 



FIG. 3: Energy, £[ mj ], of Eq. 114|l of a static pair of TSDs 
versus ln(p/a). The continuous lines are linear series ap- 
proximations: (a) when all the points are retained [£[ mj ] — 
2.15235 + 0.495788 ln(p/a)]; (b) when (p/a) < 3 are not in- 
cluded [e [md] = 2.16149 + 0.490593 ln(p/a)]. 



Table ^ summarizes the energies e[ m j ](w) from 
[m, j] = [1,0] to [7, 4] and w = 4, 6, • • • , 20. The w -> 00 
limit, £[ m ,i\i is obtained by fitting Ei m ^(w) by a power 
series in l/w. 

For moderate to long distances, our results indicate 
that the Heisenberg energy of such a static excitation 
increases as ~ hip, as predicted in Ref.l28l Nevertheless, 
a tiny deviation from this behavior is observed for small 
values of p. This is because details of the lattice are more 
important for short distances, as also is expected from the 
form of the denominator in Eq. (|14f> . See, for instance, 
Fig. |21 Therefore, it is expected that the TSDs of a pair 
will try to remain as close as possible. However, this 
is not enough to decide whether charge-wearing TSDs 
organize themselves as bound pairs of two charge- wearing 
TSDs, with bosonic character, or each charge-wearing 
TSD would bind to a spin- wearing TSD, leading to a 
fermionic-character pair. 



4>f m<j ] (k) = AW] E eik ' R E l R ' R + P[mJ])Xp [m , 

R P[m.i] 

(16) 

where Af[ mj -] is the normalization; pi m ^1 is a vector ob- 
tained by any operation of the point group acting on 
[m,j], with < j < m and m + j=odd; |R,R + p[ TOJ ]) is 
the state with a static pair of TSDs, let's say, a charge- 
wearing TSD lying on site R, and a spin-wearing TSD 
(a second charge-wearing TSD) on R + P\ m j\ for the 
fermionic (bosonic) pairs; x% 1S the appropriate char- 
acter of the irreducible representation S, with X[m,j] = 1- 
Finally, C is the square lattice (£ = A) for the fermionic 
(bosonic) pair. Then, care must be taken with the al- 
lowed values for k. For instance, when dealing with 
the fermionic pairs, k belongs to the Brillouin zone of 
a square lattice with lattice constant a. On the other 
hand, for bosonic pairs | fc^ | , | A; y | < 7r/2a, because the 
summation is restricted to run on the sublattice A. 

It can be readily seen that only the A± and B\ sym- 
metries are allowed for j=0. Since different symmetries 
do not mix, here we restrict ourselves to Ai and B\ sym- 
metries even when j ^ 0. 
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TABLE II: Energy with respect to the 
[to, j], and the extrapolation to w — > 
asymptotically equivalent to s.\ m ,j]- 
[m, j] w — 4 



half-filled ground state, E[ m j](w), of Eq. 11411 . associated with a pair of TSDs separated 
oo, £[m,i]- Values for [0,1] and [1,2] have been included to emphasize that ey, m ] is 
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14 



16 



18 



20 



[1,0 
[0,1 
[2,1 

[1,2; 

[3,0 
[3,2' 
[4 

[4,3 
[5,0' 
[5,2' 
[6 

[5,4' 
[6,3 
[7,0 
[7,2' 
[6,5 
[7,4' 



2.3044 
1.9925 
2.8537 
2.5418 
3.2688 
3.3094 
3.7245 

4.1566 
4.1636 
4.5957 



5.0307 
5.0319 



2064 
0708 
6837 
5481 
9316 
0050 
2529 
2657 
5239 
5393 
8103 

8129 
0860 
0892 



2.1730 
2.1009 
2.6325 
2.5605 
2.8071 
2.9020 
3.0766 
3.1062 
3.2728 
3.2961 
3.4923 
3.3104 
3.5003 
3.6934 
3.6992 

3.7035 



2.1587 
2.1143 
2.6133 
2.5689 
2.7508 
2.8594 
2.9969 
3.0416 
3.1525 
3.1823 
3.3379 
3.2115 
3.3521 
3.4970 
3.5056 
3.3598 
3.5159 



2.1515 
2.1213 
2.6044 
2.5742 
2.7210 
2.8385 
2.9551 
3.0115 
3.0863 
3.1210 
3.2522 
3.1629 
3.2725 
3.3852 
3.3963 
3.2897 
3.4130 



2.1473 
2.1254 
2.5996 
2.5776 
2.7034 
2.8268 
2.9307 
2.9956 
3.0461 
3.0845 
3.2000 
3.1365 
3.2255 
3.3155 
3.3287 
3.2518 
3.3513 



2.1447 
2.1280 
2.5966 
2.5710 
2.6921 
2.8196 
2.9151 
2.9863 
3.0199 
3.0611 
3.1659 
3.1208 
3.1956 
3.2691 
3.2842 
3.2299 
3.3119 



2.1429 
2.1298 
2.5947 
2.5816 
2.6845 
2.8148 
2.9045 
2.9804 
3.0018 
3.0452 
3.1424 
3.1108 
3.1756 
3.2369 
3.2534 
3.2163 
3.2853 



2.1416 
2.1310 
2.5933 
2.5827 
2.6790 
2.8114 
2.8971 
2.9764 
2.9889 
3.0338 
3.1256 
3.1041 
3.1615 
3.2135 
3.2312 
3.2075 
3.2666 



2.1334 
2.1386 
2.5792 
2.5845 
2.6509 
2.7900 
2.8617 
2.9481 
2.9304 
2.9866 
3.0553 
3.0701 
3.1078 
3.1048 
3.1366 
3.1633 
3.1980 
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FIG. 4: The hopping of a spin-wearing TSD (white star) by 
the action of Hj. 



A. Energy of the fermionic pairs 

The expectation values given by the wave functions of 
Eq. lO are 



-t (cos k x a + cos k v a) 5 



1,0] 



;[l-0] 

M 

+ J(£[m,j] +7[m,j]) 



' ,,5 



if 

kX 



[j,m 



cos k x a cos La 



where / is the ionization potential (electron affinity) 
for hole-doped (electron-doped) materials. t and t\ 
are nearest- and next-nearest-neighbor hopping integrals, 
with i = 1 when there is only one charge-wearing TSD 
in the plaquette, and i = 2 when there are two nearest- 
neighbor charge-wearing TSDs in the plaquette; 7[ m .j] 
arises from the Heisenberg terms involving the spin- 
wearing TSD (see fig-EJ), 



T[m,j] 



\x s 



3 A [0,1]' ni, — ±, 

- iX§ >TO ]^i,m-i, otherwise. 



There are two families of nonzero off-diagonal matrix 
elements of the Hamiltonian. When m — m! — ±1 and 



i-i' = ±i, 

H [Z,j] 1 = A f*i cosk x acosk y a ~\j ~ g-%m-iX§, m ] 
When m'—m — ±2 and j' — j or ml = m and j'—j = ±2, 



[m'.j'] _ 



-V3J/12, [m',f] or [m,j] = [1,0], 
-AJ/8, otherwise, 



(20) 



with A = y/2 when either j or j' is zero (but not both), 
and A = 1 otherwise. 

The zero-order lowest-lying fermionic pairs are the 
[1,0]. Close to T the energy of these fermionic pairs is 



(21) 



(17) with 



= I-t 



2t'i X[o,i] 



J 



£ [1,0] + J - nX[0,l] 



2m\, 



t t' A/ 5 



(22) 



Thence, for t[ > J/6, the zero-order lowest- lying 
band has x 2 — y 2 symmetry. This is the case for the 
Laa^Sr^CuCU (Ref. (H (LSCO). For the monolayered 
HTSC of Table 1 in Ref. 113 we know that J/6 ~ 
0.019 — 0.030 eV. On the other hand, strong differences 
on hopping integrals among the different HTSC are not 
expected, as suggested by the small variations observed 
on the nearest-neighbor hopping integral, t. Therefore, 
we expect that the zero-order lowest-lying band will show 
Bi symmetry for all of these HTSC. 
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We are now concerned whether admixing wave func- 
tions with different [m, j] to the ^ j (k) would be rele- 
vant or even if the ordering of the lowest-lying Ai and 
B\ bands could be reversed. We have obtained the cor- 
rections to Cji o] (k) and C\ia\ s ) by diagonalizing the ma- 
trix of the Hamiltonian in the basis of the two, three and 
four lowest-lying wave functions with the appropriate A 
or B symmetries. Thence, making use of the parameters 
for the LSCO^ when the number of fermionic pairs per 
Cu is p = 0.05 — 0.07 we obtain corrections to the zero- 
order Fermi energy of A^ 2 ) = —16.98 to —14.56 meV, 
A^ 3 ) = -1.45 to -1.33 meV, and = -0.26 to -0.20 
meV for the B\ band, while for the A\ the corrections 
are A< 2 ) = -7.21 to -5.1 meV, = -3.15 to -2.7 

meV, and A^ 4 ) = —0.15 to — 0.1 meV. For the mono- 
layered HTSC of Table 1 in Ref.lH assuming t\ ~ 0.2t, 
these corrections are slightly decreasing with J/t. We 
observe that, up to meV, the correction to the Ai energy 
is smaller than the correction to the B\. Therefore, it is 
expected that the lowest-lying band still has x 2 — y 2 sym- 
metry. Furthermore, the band with symmetry A\ would 
not start filling until a critical doping of p » 0.20 — 0.22 
holes per Cu02 unit, provided that all the charge- wearing 
TSDs organize as fermionic pairs. At this doping, the cor- 
rections A(") to the energies are still smaller than those 
referred above. Since this doping is out of the range of 
our interest, we restrict ourselves to consider only the 
band with x 2 — y 2 symmetry. Also, since the error in 
the parameters of the Hamiltonian are of the order of 
meV, we neglect corrections to the energy smaller than 
1 meV. Therefore, we consider only the [1,0], [2,1], and 
[3, 0] wave functions to describe the lowest-lying band of 
the fermionic pairs. 

B. Energy of the bosonic pairs 

When dealing with t-J-like model Hamiltonians, the 
electrostatic repulsion is generally neglected, although 
with some exceptions 42^i Since the screened electro- 
static repulsion between the charge-wearing TSDs in a 
pair, V] mj ], may be relevant, here it is included in the 
diagonal terms of the Hamiltonian, 

(H)f m<j] = Je [m , 3 ] +21 + V [md] + (H t + H t ,)f m<j] . (23) 

It is not difficult to obtain (H t )f n = 0, and 

( H t')fm,j] = K [™,j]Xfj, m ] (1 + cos fc^acos k y a) , (24) 
with 

[ 2t> 2 , [m,j] = [1,0], 
K [m,j] = i *ii m-j = l,m>l, (25) 
I 0, otherwise. 

The electrostatic repulsion within a bosonic pair is ex- 
pected to be screened by the gas of the fermionic pairs. 
However, the fermionic pairs have been considered so far 



as moving in a two-dimensional square lattice. It is gen- 
erally accepted that the c axis effect is simply to tune 
the electronic structure of the CuC>2 planes. Neverthe- 
less, screening is a three-dimensional effect that could be 
taken into account by an interlayer hopping integral t±. 
Considering a nonzero ij_ would imply a correction to 
the energy « t±k 2 _c 2 , where c is the lattice constant per- 
pendicular to the ab layers. Since t± is rather small^ it 
can be neglected for the energy-balance considerations, 
but it is essential for screening purposes. Then, for the 
electrostatic repulsion within a bosonic pair we take as 
a first approximation the Yukawa potential 3 ^^ as the 
dominant term, 

V p ^^-exp{~(4ne 2 g F ) 1/2 p} , (26) 

where g F is the density of states at the Fermi level of the 
fermionic pairs per unit of volume of the solid, 

/q 2 \ V3 

'*"\^wk) > (27) 

p being the number of fermionic pairs per site, and v is 
the number of square-lattice layers cutting a unit cell. 
Therefore, close to T, the diagonal terms of the Hamilto- 
nian are 

(H)f mJ] w Je [mj] +2/+^exp{-/3pM 1 / 6 } 

+ K \m,j]Xi m] ( 2 -^ ft2fc2 ) > ( 28 ) 

with 

2e / 6 \ 1/6 

The non-zero off-diagonal elements of the Hamiltonian 
can also be readily obtained, 

H [Zj] 1 = >^t'i(l + cos k x a cos ky a); \m-m% \j - j'\ = 1, 

(30) 

with A = \/2 when either j or j' is zero, and A = 1 
otherwise. 

Since the nonzero off-diagonal elements of the Hamilto- 
nian are important as compared to the differences among 
the diagonal elements, the energy of the bosonic pairs 
must be obtained by diagonalizing the matrix of the 
Hamiltonian. 

Since the screened electrostatic repulsion decays faster 
than exponentially, while £[ m ,j] increases logarithmically, 
there must be a minimum in the energy and confinement 
is expected to occur. Furthermore, from Eq. (|28|l and as- 
suming that t[, t' 2 > 0, 2 — it is expected that the bosonic 
pairs will also show x 2 — y 2 symmetry as it is gener- 
ally accepted. 3 At this point, it is worth noting that if 
the next-nearest-neighbor hopping is neglected and the 
Hamiltonian is reduced to the t-J model, it turns out 
that the A\ and B\ symmetries would be degenerate. 
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C. Two- fluids equilibrium condition 

At T = 0, the question now is whether this lowest- 
lying bosonic pairs would have lower energy than two 
fermionic pairs in its Fermi level, so the bosonic pairs 
would be favored. At low doping level, it is expected 
that the fermionic pairs will be favored. Nevertheless, 
the Fermi level, kp w y/2np/a, increases linearly with p, 
while the electrostatic repulsion among the two charge- 
wearing TSDs in the bosonic pair is exponentially re- 
duced with p 1 / 6 , so its ground-state energy is lowered. 
Therefore, we wonder whether there exist a critical value 
p c such that the ground-state of the bosonic pairs, as 
measured with respect twice the energy of a fermionic 
pair at its Fermi level, A(p), is zero. 

To explore such a possibility, we have diagonalized the 
matrix of the Hamiltonian for k = 0, and increasing val- 
ues of p. To reach corrections to the ground-state en- 
ergy within the order of meV, we have considered up to 
a 12 x 12 matrix involving all the states which would 
contribute to perturbation theory truncated to tenth or- 
der. For a certain regime of parameters, at low enough 
doping level, the fermionic pairs are favorable. As p in- 
creases there exist a critical value of doping, p c , such that 
A(p c ) = for the lowest-lying bosonic pairs. 

Doping above p c yields bosonic pairs. In this case, 
the bosonic pairs are expected also to contribute to the 
screening and the electrostatic repulsion could become 
negligible. If so, there would be a cascade process of 
pairing among the fermionic pairs until a new equilibrium 
between the two fluids is reached at pf < p c . Thence, 
we expect the number of bosonic pairs at p > pf to be 
Pb ~ (p — P/)/2. There is a lower limit of pj such that 
Ao(p/) = 0, as obtained when the electrostatic repulsion 
is completely neglected. 

For the sake of estimating the order of magnitude of 
p c and pf for a generic HTSC, let us make use of the pa- 
rameters appropriate for LSCO. We take a sw 3.8 A and 
c/am 3.47 from Ref.yl For the Hamiltonian parameters, 
we take J = 0.144 eV and t = 0.549 eV from Ref.|H and 
t' 2 = 0.130 eV and t' x = 0.112 eV from Ref. |H All of 
these parameters were obtained from high-level ab initio 
calculations with the geometry as the only external in- 
put, being the errors within the meV. Since there is no 
high-level ab initio calculation for the interlayer hopping 
integral, we use the low-doping t± w 0.7 meV obtained 
from experimental results by Zha, Cooper, and Pines4£ 
Within these parameters regime we get f3 « 7.8/a. Com- 
puting A(p) for increasing values of p we find that A(p) 
is changing its sign at p c w 0.0524 (see, for instance, 
Fig. |SJ. At Pc, the mean distance between the two holes 
of the pair is (p c ) w 9.08 A with a standard deviation 
o~ c ~ 5.08 A. Identifying in a rather loose way the spa- 
tial extent of the pair wave function (~ p c + a c ) with the 
coherence length £, we obtain £ sa 14.16 A, in good agree- 
ment with the in-plane value (£ ~ 14-15 A)2^ suggested 
from experimental findings. 

On the other hand, when the electrostatic repulsion is 



< 




0.05 0.0505 0.051 0.0515 0.052 0.0525 
P 



FIG. 5: The ground-state energy of the bosonic pairs, as mea- 
sured respect twice the energy of a fermionic pair at its Fermi 
level, A(p). Top, A(p), with screened electrostatic repulsion. 
Bottom, Ao(p), without electrostatic repulsion. 



neglected, A (p) changes its sign at pf » 0.0505, and the 
mean distance between the two holes is (pf) ~ 8.78 A 
with a standard deviation er/ ~ 4.99 A. Therefore, the 
estimated coherence length is £ 13.77 A. Again, it is 
worth noting that if the next-nearest-neighbor hopping 
is neglected and the Hamiltonian is reduced to the t-J 
model, the value of the critical doping (p c ~ 0.28) is out 
of the range where the superconductivity is observed. In 
addition, at so high doping the validity of such model 
Hamiltonians could be questioned. 

As far as we know, for other monolayered cuprate su- 
perconductors only the t and J parameters have been 
obtained from high-level ab initio calculations. Never- 
theless, for the purpose of estimating how p c and pf vary 
with J/t, let us assume that the ratios t[/t and t' 2 /t, 
as well a/3 do not change very much among them, and 
use the values appropriate for LSCO. If so, we find that 
Pf and p c increase as J/t decreases. See, for instance, 
Fig. El where we also include the pf and p c values for 
LSCO as a function of J/t. This result suggests that for 
low J/t the onset of superconductivity would be located 
at a too high level of doping such that it could be beyond 
the validity of the present approximation. Therefore, the 
superconductivity could be suppressed. 

Furthermore, it is worth noting here that the parame- 
ters that characterize a superconductor are taken as inde- 
pendent of doping. Nevertheless, J as well as the hopping 
integral do depend locally on doping, as suggested by 
the high-level emphab initio calculations of Calzado and 
Malrieui^ For instance, their calculations suggest that 
J decreases, while t increases, with doping. Therefore, 
it is expected that J/t decreases with doping, probably 
not linearly (see, for instance, the dotted lines of Fig.[f)J). 
Consequently, pb would decrease, eventually down to zero 
at a critical doping p' c such that the J/t as a function of 
the doping crosses again the pf function. Thence, the 
superconductivity would be suppressed in the overdoped 
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0.15 
J/t 



FIG. 6: Vc for the monolayered HTSC (•) and for LSCO with 
different J/t (A); pf for the monolayered HTSC (•) and for 
LSCO with different J/t (□). The continuous lines are power 
series approximations to p c and Pf sets. The dotted lines are 
hypothetical J/t vs p curves. 



regime for a doping p > p' c . Therefore, a better knowl- 
edge of the parameters is essential to fully understand 
the phenomenon of HTSC. 



V. SUMMARY AND CONCLUSIONS 

We have shown that the Heisenberg energy associ- 
ated with a pair of TSD in a spin=l/2 square lattice 
increases logarithmically with distance. Therefore, a 
charge-wearing TSD (either a hole or a doubly occu- 
pied site) in a spin=l/2 square lattice binds to another 
TSD, either to a spin- wearing TSD or to another charge- 
wearing TSD. 

We have constructed symmetry-adapted extended 
wave functions for both a fcrmionic pair of a charge- 
wearing TSD and a spin- wearing TSD, and a bosonic 
pair of two charge- wearing TSD. The energy associated 



with such fermionic and bosonic pairs has been obtained. 

For the lowest-lying fermionic band and the lowest- 
lying bosonic pairs the symmetry turns to be x 2 — y 2 
when t[ > J/6 and ti,t' 2 > 0, respectively. Since these 
conditions are fulfilled for monolayered HTSC, we obtain 
that the symmetry of the bosonic pairs is x 2 — y 2 for these 
materials, as it is generally accepted £ 

For the LSCO compound, we find a critical doping 
for bosonic pairing p c « 0.0524 and Pf ~ 0.0505 (when 
the electrostatic repulsion is completely neglected) . This 
finding could be related to the onset of High T c supercon- 
ductivity, the superconducting state being a Bose con- 
densate. This is also compatible with the existence of 
pairs above T c , a forerunner of the pseudogap physics of 
the cuprates. At the critical doping, we find a mean dis- 
tance between the two holes of the pair (p c ) « 9.08 A 
((pf) ~ 8.78 A), and an estimated coherence length 
£ « 14.16 13.77 A). These features are in good 

agreement with the experimental result of p c ~ 0.05 
(Ref. EH) and £ - 14-15 AM 

For the monolayered cuprate superconductors of Table 
1 in Ref. [12, we have obtained pf and p c as a function 
of a and J/t, while keeping fixed a(3 « 7.8 and t\/t. See, 
for instance, Fig. HjJ It can be observed that pj and p c 
increase as J/t is lowered. 

Extensions of the present work towards T ^ 0, and to 
other possible charge- wearing TSDs self-organization are 
in progress. For instance, it is of interest to explore the 
parameter regime for bosonic pairs with xy symmetry, 
phase separation, and stripe formation. 
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